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Abstract 

A sleeve valve is under development for ground-based forced response testing of air 
compression systems. This valve will be used to inject air and to impart momentum to the 
flow inside the first stage of a multi-stage compressor. The valve was designed to deliver a 
maximum mass flow of 0.22 lbm/s (0.1 kg/s) with a maximum valve throat area of 0.12 in 2 
(80 mm 2 ), a 100 psid (689 KPA) pressure difference across the valve and a 68°F, (20°C) air 
supply. It was assumed that the valve mass flow rate would be proportional to the valve 
orifice area. A static flow calibration revealed a nonlinear valve orifice area to mass flow 
relationship which limits the maximum flow rate that the valve can deliver. This nonlinearity 
was found to be caused by multiple choking points in the flow path. A simple model was used 
to explain this nonlinearity and the model was compared to the static flow calibration data. 
Only steady flow data is presented here. In this report, the static flow characteristics of a 
proportionally controlled sleeve valve are modelled and validated against experimental data. 

Introduction 

Compression systems for jet engines have long been vexed by a region of flow 
uncertainty as shown on a typical compressor map in Figure 1 . To gain an understanding of 
the fluid dynamics in this region of the compressor map, different models have been developed. 
These models need to be validated against experimental data to gain confidence in the 
predictions made by these computer codes. An experimental study of the fluid dynamics in 
a compression system requires a means by which to perturb the flow. In low speed 
compressors, (those with a rotational speed of less than about 6000 rpm) variable inlet guide 
vanes, bleed ports, the exhaust nozzle, and fuel flow have been used to excite the system. 
These actuation methods have bandwidths which are insufficient to excite some of the faster 
dynamics in high speed compression systems, (those with a rotational speed greater than about 
6000 rpm). For example it has been shown that rotating stall fluid dynamics can have 
rotational speeds of 20-50 percent of the compressor shaft speed [1]. In order to excite the 
first three modes in a circumferential modal model of rotating stall [1], a bandwidth of 3 times 
the rotating stall speed would be required. For example, a compressor with a design speed of 
18000 rpm, (300 rps), and a rotating stall at 50 percent of the shaft speed, an actuator with a 
minimum bandwidth of 450 Hz would be required to stimulate the first 3 modes represented 
by a modal model. Thus, the development of a high speed valve for injecting air into the first 
stage of a multi-stage compression system was undertaken to provide the enabling technology 
for compressor research. 
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A high speed prototype valve has been developed under the assumption that the valve 
exit air mass flow would be proportional to a variable choking orifice area within the valve. 
Static calibration of the prototype valve uncovered a nonlinear relationship between the valve 
orifice area and the exit mass flow, caused by multiple choking points in the flow. A simple 
static flow model of the valve was developed to understand the causes of this nonlinearity. In 
the following report, the static flow model is presented and the results predicted by this model 
are compared to the valve flow static calibration data. 

Valve Static Flow Model 


Figure 2 is a schematic diagram of the prototype valve, with a cutaway section to reveal 
the valve interior and the flow path. Air enters the manifold through the supply line inlet area, 
Aq. Ao is sized so that the air manifold pressure is maintain near the supply line pressure. The 
air then passes through the controlled orifice area, A,. A, is proportional to the linear position 
of a sliding sleeve that cover a set of ports in the "arm" of the "sleeve" valve, as shown 
pictorially in Figure 3. The sleeve motion has a range of ± 1 millimeter. Once the air passes 
through A„ it continues through a nearly constant area tube to the ejector exit area A 2 . The 
important item to note in Figure 2 is the relative size of flow areas A, and A 2 . The points of 
interest in the flow are: "1", which denotes the upstream supply conditions; "2", which denotes 
the lumped volume between the supply and the sink (exit); and "3", which denotes the 
downstream exit conditions. Figure 4 shows a schematic diagram of the valve model geometry 
with points "1", "2", and "3" indicated. The following assumptions are made: 


1) Aq is large enough to provide suffice mass flow to the largest orifice area. A,. 

2) The kinetic energy of the flow is lost after passing through the orifice area, Al. 

3) The system is adiabatic, (no heat exchange) and isentropic. 

4) The working fluid, air, is calorically perfect, (ideal gas with constant specific heats). 

5) No reverse flow will be encountered. 


In Figure 4, T, P, p, and w indicate the temperature, pressure, density, and mass flow rate 
respectively. The mathematical model of the valve depicted in Figure 4 is relatively simple 
and is shown in Figure 5 as a block diagram. The calculation of the mass flow through A,, 
(w,), and the mass flow through A 2 , (w 2 ), are represented by two of the blocks in Figure 5. 
The third block equalizes the mass flow through A, and A 2 using continuity and the energy 
equation. The mass flow calculations for area A, are: 


PRFj = 


w i = A i p i 
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( 1 ) 

( 2 ) 


with 

PRF, = Pressure Ratio Function (nondimensional) 
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w, = mass flow through orifice A„ (lbm/sec) 

A, = orifice area, (ft 2 ) 

P, = upstream supply pressure, (psfg) 

Tj = upstream supply temperature, (deg R) 

P 2 = pressure downstream of orifice, (psfg) 

R = ideal gas constant for air = 53.3 lb f f / (lbm deg R) 

g c = English units conversion constant = 32.17 lb f f / (lbm s 2 _) 

y = specific heat ratio for air =1.4 


Equation (1) and (2) originate from reference [2] and are derived in the Appendix. In the mass 
flow rate calculation, the pressure ratio P 2 /P, is first checked for choked flow conditions. If 
the flow is not choked then the actual pressure ratio across the orifice is used in equation (1). 
If the flow is choked, then the choked flow pressure ratio, equal to 0.5282, is used. Note that 
the mass flow calculation is idealized and no discharge coefficient is assumed. Using equations 
(1) and (2) with the appropriate substitution of "3" for "2" and "2" for "1" yields the mass flow 
rate equations for area A 2 : 


prf 2 


/p > 

r 3 

Y _ 

f p 3 ] 

ml 

Y 

P, 

L\ 2 J 


l p J 



(3) 


w„ 



(4) 


The steady mass flow through A, and A 2 is balanced by varying the conditions in the lumped 
volume, "2" and recalculating the two mass flows until equalization is obtained. This is 
achieved by using a lumped volume at point "2" and allowing continuity and the energy 
equation to provide balancing of the mass flow rates. Continuity at point "2" yields: 

2 _ ( w i - w 2> (5) 

dt " V 2 

where 

dp 2 /dt = the time rate of change of the air density in the lumped volume, (lbm/ft 3 ) 

V = the volume of the lumped system, (ft 3 ) 

An energy balance of the lumped volume yields: 

^ - YR [ T,w, - T 2 w 2 ] («) 

where T, and T 2 are total temperatures and c p is the specific heat at constant pressure. 
Equations (5) and (6) are derived in the Appendix. Since T, and T 2 are total temperatures, and 
there is no heat addition, then T 2 =T,. It can be observed from equations (5) and (6) that steady 
flow requires w,=w 2 . Eqs. (5) and (6) represent simplified volume dynamics for section "2". 
Currently, these equations are only used to balance the mass flows for this static valve model. 
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Derivation of Nonlinear Flow Relationship 


During the design of this prototype valve, it was assumed that there would be one 
choking point in the flow, that this point would be at the controlled valve orifice area, point 
"1", and that it would always be choked for sufficient supply pressures. For choked flow and 
a fixed supply pressure, PRF, from equation (1) would then be a constant, (P 2 /P 1 =0.5282 for 
choked flow). Plugging this constant PRF, into equation (2) yields an equation for the mass 
flow rate, w, that is only a function of controlled orifice area, A, and this function is linear for 
constant values of P„T„ and P 2 . As it turns out, P 2 is not constant. Under certain conditions, 
the flow can choke at point "2" and unchoke at point "1". Once the flow has unchoked at 
point "1", PRF, is no longer constant and w, becomes a function of A, and a nonlinear function 
of P 2 . The resulting nonlinear relationship can be shown in a nondimensional format by 
manipulating equations (l)-(4). During steady flow, w,=w 2 and equation (1) can be equated 
to equation (3) as shown below: 



(7) 


The upstream supply conditions, P„ T„ and downstream sink condition, P 3 , T 3 , and exit area, 
A 2 , are all known. The orifice area A, is the only variable to manipulate and the lumped 
volume conditions, P 2 , T 2 , and p 2 are the only unknowns. Noting that T 2 =T, as previously 
stated and cancelling like terms results in the following: 
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Dividing through by the fixed, nonzero values A 2 and P, yields the following nondimensional 
mass flow function: 
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This mass flow function can be used to examine the static flow trends as a function of the area 
ratio, for fixed P„ T„ P 3 , and A 2 , without regard to units. The model was run to steady flow 
state for several different controlled orifice areas, A„ for fixed values of P„ T„ P 3 , and A 2 . 
Figure 6 is a plot of the mass flow function, MFF, defined in equation (9), versus the area 
ratio, A,/A 2 , with P,=100 and P 3 =14.7 psig, T,=68 deg F, and A, =0.07 square inches). Also 
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plotted in Figure 6 are logical indications corresponding to choked flow for orifices A, and A 2 . 
For instances, in Figure 6, with an area ratio of 0.4, MFF=0.03 and both A, and A 2 are choked 
and for A,/A 2 =0.8, A 2 is choked and A, is unchoked. The choking at A 2 causes P 2 to be larger 
than P 3 . A varying P 2 and equations (1) and (3) lead to the nonlinear relationship from A, to 
mass flow. 


Comparison to Experimental Data 

The valve was calibrated in the Flow Calibration Lab at NASA Lewis Research Center. 
The mass flow was measured using a calibrated orifice upstream of the valve at several 
different supply pressures. Figure 7 is a schematic diagram of the flow calibration setup. The 
supply pressure was measured just upstream of the valve inlet area Aq. The mass flow was 
measured with the value in it full closed position to get a bypass or leakage mass flow rate. 
This leakage mass flow was used to estimate a leakage flow area in order to properly compare 
the calibration data to the model generated data which does not include a leakage model. 
When the valve is fully closed, it is assumed that P 2 =P 3 and equation (2) is used to estimate 
a leakage area, A, Ieak . A 1Jeak and w,_ leak are used as the origin for the model estimated data. 
Thus the model area ratio and the model predicted mass flow are offset by A, Ieak and w, )eak 
respectively. Figure 8 compares the actual flow calibration data to the model predicted data 
for supply pressures of 40, 60, and 80 psig, with atmospheric exit conditions and a fixed exit 
area of 0.07 inches 2 , (diameter = 0.3 inches). The model predicted mass flows do not match 
perfectly because the model is idealized, but the nonlinear trend observed in the data is 
represented by the model. It is clear that in order to maintain a single choking point at A, over 
the entire range of operation of A„ that the exit area A 2 must be properly sized relative to the 
largest controlled orifice area, A l5 if linearity is to be maintained. Figure 9 is a plot of the 
mass flow to orifice area for the prototype valve with three different exit diameters, ( 0.3, 
0.344, 0.375 inches), at one supply pressure, 100 psig. Note that the maximum possible mass 
flow increases with the increase in the exit diameter and the linear range of the valve operation 
is extended. 

A linear function for orifice area to exit mass flow makes the valve easier to implement. 
Also, the increase in useable range allows the valve to deliver more mass flow, increasing the 
effectiveness of the valve for perturbing the fluid dynamics in a compressor. While only the 
static flow properties of the valve are being discussed here, it is important to note that to 
achieve this increase in mass flow an increase in the stroke of the sleeve is required. This will 
reduce the achievable bandwidth of the valve, for a fixed orifice geometry. The orifice 
geometry can be modified to obtain this increase mass flow range without changing the sleeve 
stroke if necessary to maintain the valve design bandwidth. 
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Summary 


A static flow model is used to predict the characteristics of a prototype air valve. The 
model predicted static mass flow as a function of controlled orifice area compare favorably 
with the prototype valve static calibration data. The nonlinear relationship between controlled 
orifice area and valve mass flow is properly represented by the static flow model. The model 
can be used to properly size the valve exit area to provide for a linear relationship between the 
controlled orifice area and the valve mass flow. This linear relationship also resultsin a larger 
maximum mass flow rate capability of the valve. 
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Appendix 


Derivation of Equations QJ and (2): The following analysis of seal leakage flow was 

original done for reference [2] and is applicable here. The flow consists of two sections, a 
nozzle region and a constant area duct region. Assume that the flow is isentropic and 
adiabatic. Continuity yields: 


w = p A u 

W = T 0 p ° A M 3 


(Al) 


where 

w = mass flow rate through the nozzle, (lbm/s), 
p = density of the gas at the nozzle, (lbm/f 3 ), 
p 0 = density of the gas upstream of the nozzle, (lbm/f^m 
( ) 0 = denotes upstream conditions, 

A = duct area, (sq. ft), 
u = gas velocity, (fps), 

M = gas Mach number, (nondimensional), 
a = sqrt( yRTg c ) = speed of sound at the nozzle, (fps), 
y = specific heat ratio, (nondimensional), 

R = ideal gas constant (lb f f / lbm deg R), 

T = local gas temperature, (deg R), 

g e = English units conversion constant = 32.17 lb f f / (lbm s 2 ). 


For a calorically perfect gas (isentropic): 
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Substituting for M a" in equation (Al), and using equation (A2) to replace "T" and p/p 0 yields: 
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Using the definition of "total" enthalpy yields an equation for the square of the velocity: 
Noting that u 2 =(Ma) 2 and using equation (A4) yields: 
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Equation (A3) can be simplified using (A5) and the ideal gas law, (p 0 =P 0 /RT 0 ), yielding: 
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which is the same as equations (1) and (2) in the body of this report. 

Derivation of equation (5): This is simply the conservation of mass for a fixed control volume. 
M is the total mass in the lumped volume and V is the fixed volume. 

Derivation of equation (61: Performing an energy balance on the control volume shown in 
Figure 4, using reference [3] and assuming no heat transfer and no work done on the control 
volume results in the following equation: 
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dt 


(M) = w - - W, 


out 
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where 

cv = control volume 
cs = control surface 

e = u + V 2 /2 + gz, total energy per unit mass 
u = internal energy per mass 

Substituting for "e" and assuming the potential terms gz are small yields: 

P 2 2 

— f (u + — )pdv = - f (u + pv + — )pu dA (A8) 

dt Jcv 2 •'cs 2 


Using the definitions of total temperature and specific heats and solving the integral for the 
lumped volume and the surface integral for the two ports yields: 




(h + — ) w 
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(A8) 


Substituting for pT=P/R and using the specific heat at constant pressure results in equation (5): 
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Figure 2 Prototype Valve Schematic 
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Valve Open Valve Closed 


Figure 3 Description of a "Sleeve" Valve 



Figure 4 Schematic of Valve Model Geometry 
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Figure 5 Math Model Block Diagram 



Figure 6 Nondimensional Mass Flow Function versus Area Ratio 
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Figure 7 Schematic of Static Flow Calibration Setup 
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Figure 8 Comparison of the Model Predicted Mass Flow Rate to Experimental Data for 3 Supply Pressures 
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Figure 9 Comparison of the Mass Flow Rate versus Orifice Area for Three Different Exit Areas 


14 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMBNo . 0704-0188 

Public reporting burden lor this collection of information is estimated to average 1 hour per response, including the time lor reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of Information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington. DC 20503. 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

November 1995 Technical Memorandum 

4. TITLE AND SUBTITLE 

Static Flow Characteristics of a Mass Flow Injecting Valve 

5. FUNDING NUMBERS 
vm t cn 

6. AUTHOR(S) 

Duane Mattem and Dan Paxson 

W U— DUD— 02— DU 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E-9936 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM- 107072 

11. SUPPLEMENTARY NOTES 

Duane Mattem, NYMAInc., 2001 Aerospace Parkway, Brook Park, Ohio 44142 (work funded by NASA Contract 
NAS3-25266) and Dan Paxson, NASA Lewis Research Center. Responsible person, Dan Paxson, organization code 
2560,(216)433-8334. 

12a. DISTRIBUTlON/AVAILABILrTY STATEMENT 

Unclassified - Unlimited 
Subject Category 07 

This publication is available from the NASA Center for Aerospace Information, (301) 621-0390. 

12b. DISTRIBUTION CODE 

13. ABSTRACT (Maximum 200 words) 


A sleeve valve is under development for ground-based forced response testing of air compression systems. This valve 
will be used to inject air and to impart momentum to the flow inside the first stage of a multi-stage compressor. The 
valve was designed to deliver a maximum mass flow of 0.22 lbm/s (0.1 kg/s) with a maximum valve throat area of 
0.12 in 2 (80 mm 2 ), a 100 psid (689 KPA) pressure difference across the valve and a 68°F, (20°C) air supply. It was 
assumed that the valve mass flow rate would be proportional to the valve orifice area. Astatic flow calibration revealed 
a nonlinear valve orifice area to mass flow relationship which limits the maximum flow rate that the valve can deliver. 
This nonlinearity was found to be caused by multiple choking points in the flow path. A simple model was used to 
explain this nonlinearity and the model was compared to the static flow calibration data. Only steady flow data is 
presented here. In this report, the static flow characteristics of a proportionally controlled sleeve valve are modelled 
and validated against experimental data. 


14. SUBJECT TERMS 15. NUMBER OF PAGES 

16 

Static flow; Valve 16 _ price code 

A03 

17. SECURITY CLASSIFICATION |lB. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMrtATION OF ABSTRACT 

OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Preserved by ANSI Std. 239-18 
298-102 
























